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x/rhe re^onse of a vehicle as it moves on a rough road is 


investigated, Por constant speed motion of the vdiicle, the input is 
stationary. However some parameters of the vehicular system such as 
dashpot characteristics are nonlinear due to the inherent friction 
present, !Ehe equi-valent linearization technique is applied to deter- 
mine an equivalent linear danping vhich depend on the input to the 
system. This melhod is checked with the results obtained by numerical 
Integration of system differential equations, Por this piipose> random 
road data is generated from knovsn power spectrum of the road profile. 


I The response quantities of interest are bounce and pitch 


1 


L 




acceleiations, accelerations of seats and riders and displacaaents of 
sprusog mass relative to unsprung masses, A melhod is suggested to i 
select the parameters of the vehicular system such as to maximise the 
comfort of the rider. 


\/ ! 



CmPTER I 


INTROHJCTIOir 


1.1 GMEEiAL 

The importance of transport Tehicles and its contribution 
to modem civilization can not be overemphasized. Therefore, the 
prime objective of the automotive engineer is to provide vehicles 
TAhich are satisfactory and comfortable to the riders and to ensure 
safety of the cargo. The comfort of a rider depends on many fectors 
such as noise, smell, mental condition of the rider, ability of the 
driver etc, and vibration to which he is subjected. The cargo safety 
can be precisely defined by putting a limit on the maximum displace- 
ment and/or acceleration of the cargo at specified frequencies. 

So far no comprehensive criteria has been evolved to quan- 
titatively determine the ride comfort talcing into account the environ- 
maital conditions. Therefore, the designer has to resort to the 
evaluation of the dynamic response of the vdaicle and the rider, -vhen 
the vehicle moves on a road. The vehicle is subject to the following 
types of disturbance inputs : 

a) Internal disturbances : 

i) Engine vibration 

ii) Vibration caused due to transmission parts such as 
gear box, propeller shaft etc, 

iii) Vibration caused by rotating parts such as wheels or 
eccentricity in bearings. 
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b) External disturbances ; 

i) Input from the road 
ii) Aerodynamic forces and moments. 

Engine vibration has been ninlmized by properly balancing 
the engine. Modem mounting techniques successfiilly reduce the engine 
vibration effects on the chassis. Precision manufacture reduces 
eccentricity in bearings. There is inherent imbalance in tyres, as 
portion of the vdieel in contact with the ground, support the weight 
of the vehicle and thus gets flattened. With a tyre under suitable 
inflation pressure, this is not much of a proDlem. Under normal en- 
vironmental conditions and low speeds normally encountered on Indian 
roads, aerodynamic forces can mostly be taken care of by the proper 
design of the profile of "the body of the car. 

The major cause of vehicle vibration is due to the road 
roughness. By road roughness we do not mean the texture roughness 
of the ro3.^j but the waviness of the road. It is not possible to 
lay a perfect straight road. As the vehicle moves on these roads 
its wheels sense the waviness. The disturbing input to the vehicle 
depsids on type of road roughness as well as on vehicle ^eed. 

The excitation to the vehicle is random in nature. This 
excitation is filtered through tyres and suspension mechanian. The 
suspension chare cteristics should be arch as to produce a response 
which should not be excessive, causing fatigue and uneasiness to the 
passeioger owr b spacifiGd period of iinie# 
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1.2 IITIEATIIRE S[JRW 

Random nature of road rougliness is assumed to be stationary j 
ergodic ha-ving a Gaussian distribution. Thus it is only 

reasonable to describe the elevation of road surface statistically 
such as in terms of power spectral density. If the profiles of two 
long stretches of a road are taken and analyzed, the power spectral 
density curves in both the cases should be approximately the same to 
support the assumption of the stationarity of the profile. In prac- 
tice, the power spectral dsisity curves match at long wavelengths, 
but a large difference between the two curves is observed at short 
wave lengths. Jones has been critical of this phenomenon called 
"tail inconsistencies" [_1 3 • Srikantalah has measured runway 
roughness by two methods ; 

a) Surveying 

b) Profilometer , designed by Panchal. 

He has observed that roughness data obtained from survey- 
ing is good for long wavelengths vhere as data from profilometer 
measurement is good for short wavelengths. 

Power spectrum of road roughness when plotted on a log leg 
scale looks like a straight line ^3,4] . Macaulay has suggested to 
use two straight lines to descxribe the power spectral density. At 
higher wavelaagths the slope is higher ^3 j« Ihe discontinuity bet- 
ween two straight lines occurs at a wave laigth of ^proximately 
20 feet. This is one of the standard data used in road construction. 
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Dodds aad Robson have classified roads according to their power 
spectral characteristics and have noted the corresponding eqaations 
for power spectrum ^4 ] . 

In early fifties response of vehicle moving on a rough road 
has been investigated us iig determinastic approach [_1 } • useful- 

ness of these investigations are limited because the actual terrain 
has a random profile. The work were mainly confined to the analysis 
of simplified models operating on geometrically idealised surfaces 
vhereas the recent trend is in analysing more complex multi-degree 
of freedom systoas operating on an arbitrary terrain profile. This 
transition has been made possible mainly due to computers. 

The vduide is represented by a set of differential equa” 
tions. These equgtions can be developed either by Lagrange’s or 
Newton’s approach. The ^ston can be simulated on analog computer 
with known dynamic equations and given road conditions [_5 j • How- 
ever for nonlinear systems the pinper tool is a digital computer 
where the differential equations can be numerically integrated [1 3* 
This method is commonly loiown as time domain approach, where the 
bounce and pitch motions of the vehicle and rider are described as 
a function of time. It gives satisfactory results but great deal of 
time is spent in programming and running the simulation either on 
analog or digital computers ^6 J , 

An alternate approach is to transfer the ^stem differen- 
tial ecjuation representing the vh.icle to frequency response 
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function [_6 3 , The forcing function is then described statistically 
in the form of power spectral density. This method is Imown as 
frequency domain analysis where frequency is the independent variable. 
The differential equations are converted to algebraic equatjnns re- 
sulting in decrease of computational time. This has the added ad- 
vantage that it takes into account the randomness of road profile. 

The input to the vehicle depends on the road roughness and 
on the vehicle speed. TOien the vdaicle moves at a constant speed on 
a road having stationary property, the input to the vehicle is sta- 
tionary, Hov/ever, if the vehicle moves at a vscriable speed, the 
vehicle is subject to a nonstationary input even thou^ the road 
roughness may be stationary. Virchis [_73 J^as shown that the effect 
of nonstationarity die to commonly encountered accelerations and 
decelerations which are not more than l/2 g for road vehicles, can 
be neglected in predicting the vehicle response. The simplification 
is very much valuable because the analysis is presumably applicable 
not only to uniform accelerations and decelerations but also to 
no3aunif02?m accelerations and decelerations. 

The work done so far in the field of vehicle vibration 
assumes the system parameters to be linear. Although tyre and 
suspension spring stiffness have linear characteristic in the work- 
ing range, Ifae damping consists of linear viscous doping and 
Coulomb damping due to friction in the da^pot, Ellis suggests to 
include the effect of Coulomb damping for future work [_5 ] . In this 
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connections, ttie data supplied by M/s Hindustan Motors ltd., 

Calcutta, shows the nonlinearity of dashpot characteristics ^8] . 

Por nonlinear systems with harmonic input, one can use 
describing function approach^the -theoretical basis of which lies in 
the method of slowly varying parameters and equi-';alent lineariz£i.tion 
suggested by Krylov end Bogoliubov £9,10j . It has also been sugges- 
ted to find out equi-valent linear imping -taking into account of the 
energy absorbed in one qy-cle due to friction in a dashpot \1l} • Ibr 
random inputs these me-thods can not be directly applied. 

Ariaratanam has used ihe -theory of Markov random process 
and -the associated Pokker-Planck equation to solve nonlinear two 
degree of freedom random -vibration problems [l2] • Crandall has 
obtained exact solution of a single degree freedom system by means 
of -the Pokker-Planck equation [15] . Ihe result is con^jared wi-lh 
approximate solution ob-tained from the application of equi-valent li- 
nearization technique. They agree up to first order of nonlinearity. 
Crandell has applied perturbation me-thod -to problems in which the 
nonlinear terms in the governing equations are small £14’] . 

Puffing's type nonlinearity can be better handled in this technique. 
These methods are apparently difficult to apply to multi-degree of 
freedom systems. Even for two degree of freedom problems, in most 
cases, it is not possible to integrate the Eokker-Planck equation to 
get an exact solution. 



7 


The method of describing function has been extended to 
systems with random input Ho'l . The treatment is similar to the 
equivalent linearization technique developed by Oaughy [_15l . These 
methods belong to the class where the nature of tiae solution is 
assumed a priori. These methods are also applicable to systems with 
nonlinear parameters which are not symmetric about the origin. Equi- 
valent linearization technique like the application of EoKker-Planck 
equation or perturbation method is difficult to apply to multi- 
degree of freedom systems. 

The effect of various parameters on random vibration iso- 
lation problem can be better handled if the nimber of degrees of 
freedom can be reduced. This can be achieved by "elastic decoupling" 
and "inertia decoupling" of the system • The system with four 

degrees of freedom and two inputs is split into two subsystems each 
having two degrees of freedom with a single input. The response of 
such subsystems when subject io white noise input and ramp noise 
input have been extensively studied [^17j18 j . 

Eor multi-input, multi-output systems, the matrix of cross 
spectral densities of the response can be evaluated from the frequency 
response function and matrix of cross spectral density of the exci- 
tation [_14 ] . The response of the vehicular system should be such 
as to provide a satisfactory ride. 

Early indices proposed for evaluating ride comfort were 
incoirplete in their treatment of human response to cumulative effect 

/ 
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of fregueticies, amplitudes of accelerations, velocity, displacement 
and phase. The establidament of uniform scale for ride comfort is 
complicated because of physical nonunifo naity of test subject's 
reactions to vibration. Janeway divided the frequency band into 
three sections from 1 to 6 , 6 to 20 and 20 to 60 cycles per second{59]. 
He stated that for the lowest band the rate of change of acceleration 
or jerk -value was the criterion where as for the middle band the 
acceleration and for the upper bsnd the velocity ±L 0 Uld be -the com- 
fort criterion. In 1958 Dieckmann proposed an empirical method of 
assessing the degree of ride comfort and broke -ftie frequency band 
into three sections. He suggested acceleration, velocity and displace- 
ment were criteria for three bands respectively ^19 3 • Mitschke has 
conc]jided -that -fche acceleration ^ould be aaall for superior ride 
comfort and ride safety [193 • There are three simple criteria after 
G-oldman for subjective response to shock and vibrationi the threshold 
of perception, of unpleasan-tness aid of -tolerance [^2o3 • The latter 
two are generally difficult to identify and reproduce. 

The tests for ride comfort consists in finding the response 
of a human being to a harmonic input. But in practice "the human 
being is subject to vibrations of varying amplitudes at varying fre- 
quencies and phase. To solve such problems some guide lines are 
recommended to determine square root of the summation of -the squares 
of ride comfort indices -taken for different componaat vibrations [1 3* 
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The approach using Ihe index of ride comfort assumes that 
the contributions of vibration to ride comfort at "varying frequen- 
cies are same. This is far from truth. This is because the mecha- 
nical impedence of man -varies for different frequencies. At l 0 "w 
frequencies (l to 100 cps) the body acts as a combination of spring 
mass damper system ^ 20 ] . Por a sit -ting man the first resonance is 
between 4 to 6 cps. Between 20 to 30 cps the head exhibits mechani- 
cal resonance aad frequency range of 60 to 90 cps suggests eyeball 
resonance, Impedence depends also, on whe"ther the rider is sitting 
erect or relaxed. Similarly "the resonant frequencies are different 
for a standing man or a rider in supine position. Therefore, it is 
only natural to have weighting function as a function of frequency 
to calculate the contribution of vibration to ride comfort, 

Butlcunas translated the input spectrum into ^ectrum of 
sensation "through human body defining "the transfer function as the 
ratio of ride sensation and input spectrum ( 1 j • Butkunas's work 
is based on Goldman' s criteria of comfort. This scheme has obvious 
limitations because nei"ther "there is a single input to the rider nor 
the human system is linear. Vibration affects "the rider through 
separate inputs such as hand on steering vheel, feet on chassis 
floor and back and head on the backrest and headrest respectively. 

Ride comfort depends on the total time a person stays in 
"the moving vehicle ^21 J , Por a shorter duiation one can "tolerate 
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Tibiations of greater amplitudes. Based on this fact, the Interna- 
tional Standard Organisation (l.S.O.) approved a recoinmendation 
concerning vibration limits and human being. Por vertical vibrations 
the recommendation is sho-wn in ihe Figure (l2) vhidi indicates a 
maximum sensitivity to vibration in the frequency range of major 
body resonances. The amplitude of vertical acceleration is clearly 
frequency and time dependent, 

1,3 PRESHIT HSrVESTIG-ATION 

It has been emphasized in earlier publications [_1>5 j that 
the effect of friction in dashpots, used in the suspension mechanism 
of a vehicle should be taken into accouat. Prom the experimental data 
supplied by M/s Hindustan Motors Ltd. f 8 ] , it is evident that fric- 
tion gives rise to nonlinear characteristics of the dadipot response. 
Thus the differential equations governing the dynamics of the vehicu- 
lar system becomes nonlinear and coupled. Equivalent linearization 
technique using the minimum mean' squared error.. criteria is applied 
to determine ihe equivalent linear parameters of the system so that 
the system will give approximately same response as would be given by 
a nonlinear system. 

The input to the ^stem is described by a power spectrum. 

As data for Indian roads were not available, the data of a runway 
is used to calculate the power spectrum. Dodds and Hobson £ 4 j have 
described different types of roads based on their roughness and their 
corresponding power spectrum. These are used as the input to the 
system. 
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A model for the -vehicle is developed. The model has 
several parameters. The systea is decoupled statically and 
cally to two subsystems. Frequency response functions are determined 
for 2 and 6 degrees of freedom systems. With kno-vm. complex frequency 
response and input power spectrum, the power spectrum of different 
response quantities is determined. The response of the seats in a 
vehicle act as input to the rider. 

Human being act as a linear system under low frequency vibra- 
tion conditions {^20^ • The discomfort of human being is frequency 
dependent. The resonant frequencies depend on posture of the rider 
such as sitting, standing etc, A transfer function of human system 
has been derived from the I.S.O, recommendation for effects of vibra- 
tion on man. Using -this frequency dependent transfer function the 
response of the rider to bounce and pitch motions of the vehicle are 
determined. The response gives an idea of whether the rider is on 
threshold of perception or in the intolerance range. Comfort is a 
function of the r.m.s, acceleration of the subject. 

There are several system parameters v\hich can be varied to 
obtain an optimum condition resulting in best possible comfort to the 
rider in a vehicle moving over a given rough road. The criteria for 
selecting Ihe system parameters is to have minimum r.m.s, acceleia- 
tion of vibration of the rider subject to a limited wheel excursion. 

The results ob-tain ed by using equivalent linearization 
technique is verified by numerical integration of system differential 
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equations, For this purpose random data is generated from a knom 
power spectrum. These data are used as inputs to the system. 

Gill's modification of Runge Kutta method is used for numerical in- 
tegration purpose. 

The thesis ends with the discussion of tae problem. 

Several conclusions are drawn from this investigation. Some sugges- 
tions are made for future work. 



OH&PTER II 


ROAD ROUGHNESS 

2. 1 SPECTRAL DECOMPOSITION OP ROAD ROUGHNESS 
2.1.1 Intro duction 


The autocorrelation function R (t, , t„) of a random 

XX 1 d. 

process X (t) is defined as 

(t^ , tg) = E [X (t^) X (t^) ] (2.1) 

Por a stationary random process the autocorrelation func~ 
tion is a function of difference of time t^ and t^. By definitionj 
this function is symmetric. 


R 


XX 


(r) 


E 


XX 


( - r) 


( 2 . 2 ) 


Ihere X = '*'1* ^ arbitraiy function h (t) the autocorre- 

r '1 

lation function satisfies j 14 j . 

A function such as R (t^j t_) vh-ich satisfies equation (2,3) is 
said to be nonnegative definite. A theorem due to S. Bochner ^14 1 
asserts that every nonn^ative definite function has a nonnegative 
Pourier transform, if such a transform exists, that is, 


S (C.3 ) 

XX 


2 IT 




R iZ) exp (~ luy'C) dr ^ 0 

XX 




(2.4) 


Prom the Pourier inversion formula it follows that R ( ’?^ ) can be 


XX 


expressed in terms of S ( O J ) : 

jCa. 

f'CO 




S (CP) exp (iC*-‘''r) du) 

XX 


(2.5) 
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Equations (2,4) and (2,5) are popularly Inom as Weiner-Khintciiiae 
relations, A physical meaning can be givoi to S ( cC- ) by consi- 
dering the limiting case of equation (2,4) in which X. = 

(0) = E (_ = j S (cc) d; :• (2.6) 

The mean square of the process ecpals to ihe sum over all frequen- 
cies of S ( ) so that S ( (O) can be Interpreted as a mean 

XX XX 

square spectral density. Since the electrical engineers Yvrere the 
first to usefully exploit this, ihey iateipreted S (<- 0 ) as po?/er 
contained in a wave and gave an alternate lame as power spectrum. 

This concept can be analogously carried over to mechanical systons. 

L 2 

Z (t) j is 

proportional to the average strain energy and if Z (t) represents the 

velocity, then E Z^ ("t) J is proportional to the average kinetic ; 

i 

energy and S (hJ ) describe the distribution of the total mean squ- i 

XX I 

I 

are "value over the frequency domain, ' 

Spectral density S ( (a? ) is an even function. In prac- | 

XX ^ 

tice frequency is only positive and esipressed in cycles per second ' 

I 

rather than radians per second. The esiperimeatal spectral density is j 

1 

denoted by G- (f). The relation between (f) and S ( OTJ ) is i 

XX XX XX [ 

t 

given by, j 


G 


it) = 4 tr ( ^ ) 


In place of (2.6) we will have ^ 
(0) - E [ ] = 


,00 


G_ (f) df 

XX ' ' 


(2.7) 

(2.8) 


/ 

o 
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2,1.2 Measuraaent of Road Roughness 

Prior to deciding iiie meliiod to he used fer measuring the 
road roughnessj it is necessary to imov/ the frecjiency range of in- 
terest for power spectrum calculations, Por the human vibration 
problem the frequency range of interest is 1 cps to 80 cps. If the 
Thiele -vibration problem is investigated for the vehicle moving at 
a speed of 30 to 75 mph. , the wavelergth of interest can be calcula- 
ted from the following relation : 

f = V / ^ (2»9) 

vhere v = vdaicle speed (ft,/sec, ) 
and ^ = wave length (ft.) 

Iter the frequency range and v^icle speed mentioned, the 
minimum wave length of interest is about l/2 ft. and maximum 110 ft. 

One of the most commonly used me-thods of measuring road 
surface is surveying by level and s-fcaffs. If sampling, which is de- 
fined as the points at which data are observed, is too close together, 
it will yield correlated aid hi^ly redundant da-fea. This results 
in wastage of time in sur-veying as well as computation. If sampling 
is too far apart, -the points will con -tain the ampliixides not only of 
the wave lengths, we think we are analysing, but of other harmonics 
as well. This phenomenon is known as ’’alising". If the roughness 
in "the higher harmonics is relatively anall, this need not matter 
gr^tly, but if the roughness is relati-vely large it can cause 
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misleading results. A sapling theorem credited to Shannon states 
that the hipest frequency that may he observed hy sampling a wave 
at a discrete inter'val A is equal to l/2 , known as lyquist 

frequency ["1, 22 [ . Therefore^the Interva.! that should b e employed 
for sampling is equal to l/4 ft in order to avoid alising. 

The sampling interval of l/4 ft. is too close to be obtain- 
ed in surveying. Therefore^ for most of the roads and airfields run- 
way roughness survey up to date, the interval is 2 ft. |" 22 j , 
Obviously this method gives an incomplete picture of profile rough- 
ness and has to be suplemented by some other method. Many methods of 
meaairing road roughness have been listed by Macaulay | 3 | • Srikan- 
taiah f 2 ] has used Panchals profilometer as a suplement. 

2,1.3 Oalculation of Power Spectral Daisity lUnction 

The random road profile rou^naess, besides being stationary, 
is assumed to be weakly ergodic. This implies ths.t temporal averages 
can replace the ensemble averages. This is very useftil because from 
a single long record of roughness data, the mean and autocorrelation 
function can be calculated. This avoids the need for collection of 
a huge set of data, Por calculation of power spectral density func- 
tion, the autocorrelation flinction need to be computed first, 

Por N data values i Y * ! , n = 1, 2, ..., N, the 

I 33. ’ 

autocorrelation fhnction, E , at the displaconent rh is calculated 

[ 23 1 = 



_J 

H - r 
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R 


r 


R (rh) = 

y ^ 


R-r 

Y’ Y' 
n n+r 


r— Oy 1y 2y • 0 «y 

(2.10) 


Where r is lag number, m is marimum lag number, h is the sampling 
interval and Y' is the road height at a point n , measured OTer a 
certain mean. Y' is calculated from the follov/ing relation. 


Y’ = Y 
n n 


J. 

N 


n’ 


n = 1, 2, 


• . • , 


R 


(2.11) 


m 


Where Y^ is the road height read off over a datum at a point n . 
The second term on right hand side of equation (2.II) denotes the 
mean of j Y^ j , n = 1, 2, ... , R. 

A procedure recommended in spectral analysis technique is 
prewhitening ^ 23 | . The process of prewhitening amounts to applying 
a fecial digital filter to the data which wall result in the fil- 
tered data having a flat spectrum. The digital filter used is given 
by, 


(2:) = Y^ (X) - (Z - h) (2.12) 

The initial or "raw" estimate of the power spectral density function 
are calculated from the Allowing relation ; 



for k = 0, 1, 2, m (2*13) 

Where 0^^ = 0 ( ,0. ) } O. = , k = 0, 1, 2, ..., 

li = corresponds to cut off frequency. 


m and 
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The raw eatimate of ihe power spectrum are then smoothened by 
’’hamming'* |* 24 "j . 

0 ^ = .54 0 ^ + .46 0 ^ 


0 . 


.23 + .54 4 + .23 


0 = .54 + .46 f( , (2.14) 

m ^m '^m-l 

A check sum is con5)uted to check the computations of the esticetes 

M •• 


CHKaJM = 


■n 


mh 




, m-1 

+ 0 ) + y 




(2,15) 


To compaisate for preiAhitening, the smoothened spectium is recoloured 
by -tie following operation ■ J 24 "| : 

0 ^ = 0 ^/^ (l-cosD.h) (2.16) 


2,1.4 Eesults 

The power spectral density is not plotted against wave 
length but spatial frequaicy 0- (radians/ft) . This is useful be- 
cause to obtain the input to the vduicle in radians per second, it is 
only necessary to multiply the road frequency by the vehicle speed v 
(ft./sec.). Similarly the amplitude of power spectral densiiy func- 
tion is expressed in ft /rad/ft., since to obtain the aveiage ampli- 
tude in any frequency band of interest, it is only necessary to 
measure the area under the curve between the upper and lower frequency 


limits 
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!I!he data for Indian roads are not available. To calculate 
the power spectral density function, data is taken for a taxi way 
£ 22 ] . The data is given for 1501 points at an interval of 2 ft. 

The power spectral density is calculated and plotted in figure ( 1 ). 

The power spectrum is also calculated for ihe two 1500 ft. inter- 
vals and also plotted in figure ( 1 ). 

2.1.5 Discussion 

It is seen f rom figure ( 1 ) tliat the power spectrums of the 
tvro stretches of road profiles are quite close at low frequency range. 
The difference becomes large at higher frequencies. This is descri- 
bed by Jones as- "tail inconsistency" fl ] • fact that power 

spectrums are almost same in the large range of frequencies confirms 
stationarity of the road profile, from -fche plot of power spectrum 
on a log-log scale (figure 2), it is observed ttiat the power spec- 
trum is approximately a strai^t line. 

Since the measurements available are at 2 ft. interval , 
consistent result is expected upto a -wavelength of 4 ft. (using 
Shannon's Theorem). But for the s-fcudy of vehicle vibration, we need 
the measurements of the road profile at l/4 ft. intervals as calcu- 
lated earlier in article 2.1.2. Thus the power spectrum plotted in 
figure ( 2 ) is incomplete, so far as study of vehicle -vibration is 
concerned and the da-ta cannot be used for the investigation. If a 
large set of road profile data is a-vailable using a calibrated 



owcr Spectral Cen 
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profilometer, it can be effectively used for calculation of vehicle 
response. 

Profile data for Indian roads are not available. Therefore 
the power spectrum data for this investigation is taken from the one 
given by Dodds and Robson T 4 1. They have surveyed different European 
roads and have classified them according to roughness. The power 
spectrum consists of two strai^t lines meeting at a wavelength of 
6,3 metres and are defined by the formula, 


! 

1 

j 

1 ® 

-B 

(n / n^) 

n < n 

0 



i ® (“o) 

(n / 

'O' 0 

(2.17) 


vdiere = 1 / 6,3 (cycles/m, ) 

2 

G (n^) is the roughness coefficient (m /cycle/m) 

Bg are slopes of power spectrum on a log log plot. 
The power spectrum data for different roads are given below 


Road Class G(n ) / 10”^ ^2 

Mean S.D. Mean S.D 


Motor ways 

Very Good 

2-8 


Good 

8-32 


Very Good 

2-8 

Principal 

Good 

8-32 

Roads 

Average 

32-128 


Poor 

128-512 

Minor 

Average 

32-128 

Roads 

Poor 

128-512 


Very Poor 

512-2048 


1.945 .464 1.360 .221 


2.05 .487 1.44 .266 


2.28 


534 1.428 


263 


2*2 CORREIA TIOH OE THE INPUTS 


The vdiicle has inputs depending on the number ox #ieels. 
By the assumption that both front and rear Vv^heels are subject to 
the same road roughness^ tlie number of inputs becomes two for a Te- 
hicle with four wheels. If it is further assumed that rear wheel 
follows the front wheels then the input spectrum to both the viieels 
ax'e the same and at a particular instant inputs to the two wheels 
are related. The relation depends on the wheel base and vdiicle 
speed. Clearly the input to the rear wheel lags the input to the 
front wheel ty a time given by {6 ] , 


Ik T 


( 2 . 18 ) 


= 0 i = k 

Where d is wheel base (ft.) 

The cro ss— correlation function for the two inputs arej 
By -the application of (2.4) we get the cross power spectral density. 


(2.19) 


S (<>■'') = 

X JT ^ ' 


.-riO 


> , . “1 


^ ( t ) e3cp (-iOU^C) d'Z 


i>0 






letting t = ‘T - 


5k 

- j (t) exp (-±cM) exp (~i 6J ^j^.) di 
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Substitutnng equation (2.5) in (2.20) 




(2.21) 


2.3 EimTIOFSElP BETVTEH^ DiSPIACaiMIi: , l/EIOCTiy Mp 
AGCEIERATION POYffiR SPB 3 TRU 11 JS OP THE INPUT 

Eor a stationary random process, 

^ ) = R^ (t - s) = E [ X (t) Z (s) J (2.22) 

Di£ferentaatiag ttie eq.uation (2.22) with respect to t and s and 
noting that differentiation and 635)6 ctation commute, 

a-r^ hx = aF^ ® t 

= B { X (t) i (b)] = By ( ^ ) < 2 - 23 ) 

Differentiating equation (2,5) with respect to t and s 


2 2 

^ ^ a^(o3)expfir.= (t-s)] a<u 


-'-GO 

.2 2 

— **■ i Qfj I S„ 
P ^-00 ^ 

~ (.0 ^ ^ ^ 


^XX exp I i CO (t - s)jdu 

: ) (2.24) 


As in equations (2.4) ariu (2.5) the Weiner-Khintchine relation for X 


can he written as 


S” (Gu) = h^x (-icor) dT 


(T) = I (lO ) exp) (ico't) dco 


(2.25) 


By taking inTerse Iburier transform of equation (2.23) and using 
relations (2,24) and (2,25), 
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(Ck' ) exp (i oC' r) 


(tja) = 


Similarly, it can. be shovm itLat, 


^ Sjj^ (o5) exp (ic^'’") doJ 


(2.26) 




(2.27) 


2.4 GHTEeAIIOfr OP EiNlDM DATA 

To generate random road data, at some specified inter'vals, 
■whiob coixesponds to a given power spectrum, we use a method as 
described below f 25 1 , 

V . s 

Let 1 (t) = ^ A^ Cos (2 IT t + 2 rr 4) (2.28) 

k=1 

Hiere 0, is the Kth sample value of uniformly distributed random 

ilC 

■variable in 'the inter-val (o , 1 J that is, 


f0 (x) = 1 if X ^ (0, 1] 

= 0 , otherwise 

fjj. = (k - ^ ) A f 

and are constants depending on 'the power spectrum. 

N i~ r 

E [ X (t) "] = A^ E Cos I 2 11 (^ "^ + J 

k=1 


(2.29) 

(2.30) 


[ Oo 3 ^ 2 TT j ] 

1 'I 

= f Cos I 2 TT (:^ t + I d ^ 
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Therefore, E | X (t) 


I X (t) X (t + C ) 


= 0 (2.31) 

IT 

= ^ ^ A. A E i Cos ; 2 fl (t + T) 

J=1 k=l ^ ' 


+ 2 TT ( Cos f 2 TT t f 2 TF j2f. J | 


How let. 


I = B j Cos J 2 n f^^ (t + £; ) + 2 ]1 I Cob j 2 v ^ + 2 jr 0^J ~j 


Case I : 


If j = kj 


^ 1 

I = J Cos ^ 2 Tf f . (t + t ) + 2 If 0^ { Cos I 2 17 f^ t + 2 TT d 0 


(2.32) 


Case II : 


If 3 0 k, 1 = 0 


Therefore, E ^ ("t) X (t + X )J *= Ej; ( "^ ) 

H A.^ Cos (2 IT f Z ) 

„ ^ -J 3 

- 2 
3=1 

Weiner - Haintclaine relation^ip can be written in the form, 


(2.33) 


(f ) = 4 ( t ) Cos (2 IT f T ) d "C for 0 < f < 

^ ^ (2.34) 

V (T) = f G-. (fl Cos (2TTf T) df (2.35) 


) = j G^ (f) Cos (2 IT f ‘T ) df 

Irom equations (2,33) and (2.35), 

E A ^ 

^ -_i_ Cos (2 T f . T ) = (f) Cos (2 TT f T ) df 

d=i 


(2.36) 
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Or 

Or 



= (f^ A £ 

= [2 (f^) A f ] 


(2.37) 


Af must be small, such that A. will represent the power spectrum 

D 

correctly, let be the cut off frequency beyond which the one 
sided power spectrum is insignificant. U can be found cut from ihe 
following relation ; 

N A f = fg (2.38) 

Using equations (2.28), (2.29), (2.30) and (2.3?) the random road 
data can be generated for a loiown power spectium. 


OHAPTEE III 


MA.THEMATEAL MODEL OE THE EEHIOIS 

3.1 INTEOnJCTIOH 

An Idealized mathematical model of the vehicle is necessary 
ia order to study the dynamic response of the vdiicle subject to 
random road input. Representative model has obvious advantages over 
the alternative approach of manufacturing the vehicle and then con- 
ducting physical tests, A model vhich is too simple can hardly pre- 
dict the actual response of the system. On the other hand a highly 
complicated model becomes very difficult to analyse. Ihe automotive 
engineer has to draw a compromise and select a model representing 
the vehicle for the determination of responses, 

3.2 ASSUMPTIONS 

Certain assumptions are necessary to make a mathematical 
model from the physical system. The assumptions make the mathematical 
model easier to handle and approximately represent the physical sys- 
tem in the range of inputs that are usually expected in the systan. 
These assumptions are ; 

a) The road is hard and does not deform when vdiicle 
passes over it. 

b) There is no loss in tyre contact with the road. If the 
vehicle loses contact with the ground, it :fells to the 
ground due to the gravitational force. Eor loads, en- 
countered, loss of tyre contact with •he road is not a 


frequaat phenomenon. Under iiiis assumption the vehicle 
tyre follows the road profile closely and gravitatioiBl 
force does not play part in the governing differential 
equations. 

There is bounce and p itch motions of the vehicle due to 
different inputs at front and rear vheels. Even for ihe 
front two wheels the input may be different. But the 
difference is small compared to the e2q)ected difference 
of inputs between -the front and rear wheels. Another con- 
tributing factor vshich suggests to neglect the roll motion 
is the fact, that the suspension springs and da^pots are 
not located directly above the wheels, but tovvards the 
centre of the vehicle. 

The chassis is assumed to be rigid. 

The damping in tyre is neglected. 

The tyre and suspaision stiffness character sties are mea- 
sured by static tests. The stiffness is assumed to remain 
the same under the dynamic load conditions, Ellis [^5] 
has observed that the dynamic spring constants are less 
than the static spring constants, which implies that the 
natural frequencies at dynamic condition is less than 
that at the static condition. 

The pitch axis is assumed to remain stationary with res- 
pect to the chassis. The vehicle has the pitch motion 
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3.3 BASIC EQUATIONS 


The inputs to the model as daowi m figure (3)} a3?e 
vertical road displacements and TBe bounce and pitch motions 

are affected by the load (number of passengersj baggage etc.) on the 
vehicle. Horizontal disturbing forces due to wind gust etc. have 
negligible effect on bounce and pitch motions. The relationship 
between vertical displacements ^ Y^ of the support points and 
bounce and pitch motions are given by, 

Y, = Y + I. 6 ^5.l) 

1 1 

Yg = Y + Lg e (3.2) 


Yheve - a 
and Ijg = -’b 


The resultant forces and moments due to forces on the supports j 

^2 

■Jg - I'l + ^2 ^2 ^ 

Eesolving the forces acting on the sprung mass and talcing moments 

about the centre of gravity, 


MY = f ^ - f, 
J © = <la " % 


(3.6) 

(3.7) 


Eor the front suspension, 


M 


1 






vtiere, 

= C^ (X^ - E^) 
and f!j = (^1 - 


(3.8) 


( 3 : 9 ) 

( 3 . 10 ) 



A/VW 


I 



FIG 3 MODFl OF THE VEHICLE 



Similarly, for the rear suspension, 


* 2 % - ^ 2-^2 
’siiisre, 

fg = (Xg - Eg) (3.12) 

and = Sg (Eg - Yg) + Eg ( Eg - Yg) (3.13) 


The equations of motion for the seats and passengers are given by^ 

(t) =-K^ (Z^ - Y - 1^ &) - - Y - 1^ ©) 

(3.14) 

and mg *Zg (t) = -Kg < Zg - Y - Ig &) - dg ( Zg - Y - Ig 4) 

(3:15) 


The force and momaats due to loads are givaa hy, 

f^ = - K^ (Z^ - Y - 1^ ©) - d^ ( Z^ - Y - 1^ e ) 

- Kg (Zg - Y - Ig &) - dg ( Z^ - Y - Ig 4 ) (3.16) 

and a-j, = - (Z^ - Y - 1 ^ 9 ) - 1 ^ d^ ( - Y - Ig &) 

- Ig Kg (Zg - Y - Ig e) - Ig dg ( Zg - Y - Ig 9 ) 

(3.17) 


Substituting equations (3.9) and (3.I0) in (3.8), the equation 
of motion for the front unsprung mass is obtained as ; 

(X^ - EP + ( Y + 9 - Ep + (Y + 9 - 

(3.18) 


Similarly, for the rear unspruiqg mass, 

Mg ^2 = Cg (Xg - Eg) + Sg (Y + Eg 9 - Eg) + Eg ( Y + Eg 9 - R 

(3.19) 


The equations for pitch and bounce motions are similarly derived 
They are : 
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MY = - (y + e - ( Y + ) 

- s^ (y + I2 © - R^) - ©2 ( Y + Lg © - Rg ) 

- (y + 1^ © - Z^) - d^ ( Y + 1^ © - ) 

- Kg (Y + Ig © - Zg) - dg ( Y + Ig © - Zg ) (3.20) 

and J © = - s^ (y + © - R^) - 1^ ( Y + © - R^ ) 

- Sg Ig (Y + Ig © - Rg) - Dg ^2 ( ^ + ^2 ^ “ ®2 ^ 

- K^ 1^ (Y + 1^ © - Z^) - d^ 1^ ( Y + 1^ © - Z^ ) 

- Kg Ig (Y + Ig © - Zg) - dg Ig ( Y + Ig 6 - Zg ) (3.21) 


3.4 DECOUPIIHG THE SYSTEM REEHIREtTTlAI EOJATIOKS 

The governing differential equataons derived in previous 
article consist of a finite number of system parameters which can 
be varied. If the number of degrees of freedom and number of system 
parameters can be reduced, the vibration isolation problem can be 
handled easily. This can be done by decoupling the system. 


laplace transfoim is applied to "the basic differential 
equations (with aero initial conditions) and rearranged. A block 
diagram representation for the system is shown in figure (4) -where. 



(3.22) 

(3.23) 


Sg + s Dg 
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Since not all the matrices are diagonal, the noninteraction condition 
can he satisfied if the product of the three matrices £ h^”| 


is diagonalized. That is : 

. I- 

i_ + - 1 - 
2 _ 2 


1 


M s 

1 




J s 


l 1 12 


M s 

1 


J s^ 


2 ■ 2 
s J s 


Ms J s 

Clearly, the condition for diagonalizaticn is 

1 ^1 ^2 

-JL^ == 0 

Ms J s 

Substituting equation (3 • 3 ) (3.20, 


should he diagonal 


(3.24) 


J = M ah (3.25) 

o 

Since J = M K , the front and rear suspensions will he Independent 
with regards to their inputs and outputs if, 


R = / ah = 1 (3.26) 

This condition is Inown as "inertia decoupling", -^idaich 
permits the distributed mass M to he replaced hy two point masses 
situated Tertically above the axles. The two support points become 
reciprocal oentresof percussion# 

If E (s) = R (s) = 0 , then from the block diagram 9 

^ i 2 



Since the matrices and I are diagonal, the condition 

for noninteraction of bounce and pitch motions of the body with 
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respect to the force and momeat is that the product 

of the three matrices is given hy s 


H. 


H 


, i 

r ' 


is diagonal. By actual multiplication the product 




*^2 ^2 




In order that this matrix he diagonal ; 

G I + G„ I =0 (3*28) 

11 2 2 

Substituting equation (3«3) tn (3,28), 

a G, = b G (3.29) 

1 2 

Equation (3.29) iii 5 )lies, 


^ i. = 1 (3.30) 

®2 ®2 ® 

Under these conditions, the force f^ excites only the 
bounce mode viiiile the couple q^ excites only the pitch mode. This 
condition also ensures equel static deflection of the suspension 
at the front and the rear. The system is said to be completely 
noninteracting if equations (3.26) and (5.30) are satisfied simul- 
taneously. 


3.5 EREQUEHOY EESBOISE lUNCTION EOR TWO DEGREE 
OF FREEDOM SYSTEMS 

Tn the case vhere two masses are connected as in figure (5) 
to an oscillating foundation by linear springs and linear dashpots, 
the equations of motion in terms of absolute displacoaents are : 



FIG'S TWO D'FG^’FF OF r"'-£f:r;CN" 



F"IG. 6 DAMPING FORCE VELOCITY ACROSS 

■ nonlinear DASHPOT 
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m w = -C(w-Z)-s(W-Z)-D(W-Z) (3.31) 

O 

and Z = s (w - Z) + D ( W - Z ) (3.32) 

In order to standardize the siihsequent treatment, the following 
notations are introduced 

= y G / J B / I^r 

f = /m^ 52 = ""Z " 

The relative displacaaents between #ieel and tcriaia (RW) and chassis 
and axle (RS) are, 


EW = W - X 
RS = Z •- W 

Substituting the equations (3-33) and (3.34) ia the equations 


(3.34) 

(3.31) 


and ( 3 . 32 ), 

EW + X* + EW - / 2^ ^ ^ -^2 ^2 ~ ° 

and EW + RS + X*+ RS + 2 12*^2 = 0 


(3.35) 

(3.36) 


There are four different responses narcBly the two relative 
displacaaents defined in equation (3.34) and two absolute displacaaents 
of the two masses. The fixst pair is obtained by setting [l?] , 


iwt 

X = e 

TT 

EU = 

i Wt 

and RS = ® 

in equations (3.35) and (3.36). 


(3.37) 

(3.38) 

(3.39) 
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inhere and are "Sie complex frequency responses. 


Solving these, 


Where 


( ul? ) _ - / /\ 

A = - 2 i (l + f ) '52 ^'“^'2 

- CO^ + (1 +f ) dJg I 2 i 32 

Substituting equations (3.40) and (3.41) 


(3.40) 


(J) 0-^2 OJ^ 

(3.41) 

in (3.32), the 
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frequency response function for the absolute acceleration Y can be 
obtained as ; 


= ( 2 i ^2 60\ ) / A (3.42) 

The absolute acceleration of the chassis and the displacement of 
chassis relative to the axles are of interest. Their frequency res- 
ponse functions are given by equations (3.40) and (3.42), The fre- 
quency response functions are useful in obtaining the power spectrum 
of the output if the input power spectrum is kiown. 


3.6 EREOIMCY RESPCMSE OTCTION EOE IlHEA-R STSUCTUEE 
ITETH PIEITELY MMl DE&REES OE TREBIOM 

} 

In the case of one or two degrees of freedom systems, the 
frequency response function between inputs such as displacement or 
acceleration and outputs such as absolute acceleration or relative 
displacement are simple and direct. But in ihe- case of finitely 
many degrees of freedom systems it is not as simple. In dealing with 
the latter case it is convenient to use matrix notation. Therefore 
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the dynamic characteristics of a linear finitely many degrees of 
freedom can be specified by a ma-brix of frequency response function 
[h (40)] . 

To explain the physical meaning of a typical element 
(CO) in the matrix j^H ( 'uO ) j , the follovdng is considered. Let 
a single sinusodial excitation act upon the system at a given point k 
and in a specified direction so that excitation is described as 

(t) = A exp (iCOt), -vifoere A is a complex coefficient. It is 

understood that the real part of tie complex quantity only represents 
the physical input. The steady state response at point j has same 
frequency as the excitation (t) = A ( <-^ ) exp (i (Dt), 

Therefore, the frequency response function ) is the ratio of 

steady state response at point j to the sinusodial excitation at 
point k. In general this function is complex. For physical systems 
the elements in the matrix of frequency response function satisfy 

The equations (3.14), (3.15) and (3.16) to (3.2l) can be 

represented by a single matrix equation s 

[m] X + [c] X + [k] X = f 

■where X = ^2 ^ ^ ^1 ^2 j 

Since 6^ = ’ 


(3.43) 
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0 0 0 0 0 
0 0 0 0 0 



0 M 0 0 0 

0 0 J 0 0 


0 0 0 0 0 
0 0 0 0 0 






f 



0 

-^1 

-L D, 

1 1 

0 

0 ' 

0 

^2 

-^2 

— T 1 T) 

2 2 

0 

0 

-^1 


D^+D2+d^+d2 

1 ^ D ^ ”^^' 2^2 1^1 ^^2^2 

-^1 

“^2 


"^2^2 

^1^1-^W^1^^l2d2 


“^ 1^1 


0 

0 

-‘^1 

-l^di 


0 

0 

— ■ 

0 

CM 

1 

“^2^2 

0 


0 +S. 

1 1 

0 

"®1 


0 

0 

0 

C 2 +S 2 

-®2 

“^2^2 

0 

0 


-Sg 

s^+S2+ki‘tk2 

^1^l'^V2‘^^1^1 

+ 1215:2 

“^1 

-1^2 

1 1 

"^2"2 

V 2 

+I^k^+l2^2 

+l^^+l|k:2 

-^ 1^1 

*^2^2 

0 

0 

-^1 

-kll, 

Ic ^ 

0 

0 

J 

0 

"^2 


0 

^2 ^ 


( 

= ) 2 ^ ° ° 


0 0 
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The matrix of frequency response functions , j^H ( ch )jcan 
be obtained by letting the excitation vector f = f ^ exp ( ioJt), 
\There f is a vector of constants, let the response vector be, 


H 


r 


i H ( u") ) 


( ) I f . Thus j 

= fml + i Oil C ] + j k 1 I 


-1 


I ^ 

Since j m J , [c j and [k | are 6x6 matrices, ( tO ) j is also 
6x6 matrix. 


( 3 . 44 ) 


a 


3,7 METHOD OP EQJIVAIENT LINEARIZATION 


The dBinpiag coefficients and spring constants considered in 
section ( 3 . 3 ) and ( 3 . 6 ) are nonlinear as ia shown in figures (?) io 
( 11 ). The method of equi-valent linearization has bean used as follovss 

The differential equation of motion, for nonlinear oscilla- 
tor subject to stationary Gaussian excitation, is given by , 

X + g ( X , i, t) = f (t) ^3.4 

■vvhere g ( X, X, t) is a nonlinear flinction. Further g ( X, X, t) 
is assumed to be weakly nonlinear ij5 j • 


Eowriting equatiaa (3.45)} 

X + S X + Ijr + e (^» ^ 

vfaere D is the''e(Siivalent linear damping” coefficient per unit 
mass, S is the "equivalent linear stiffness" coefficient per unit 
mass and e (X, t) is the error due to linearization. 











Mean Vei. Ins /Sec. 

' DAMPING FORCE CURVE FOR HYDRAUUC5 MAKE 
/A (AT 7-1442) (Ccurtsey : Hindustan Motors Ltd. Calcutta)' 
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If e (z, Zj t) is neglected, equation (3»46) becomes linear 
and may be readily solved. The logical choice of S and H is 
therefore, those "values that make e (z, Z, t) a minimum. One vvay of 
achieving this is to mininu ze the expected mean— squared-error . 

Comparing equations ( 3 . 45 ) and (3.46), 

e (Z, Z, t) = g (Z, X, t) - D Z - ( 3 . 47 ) 

Uow minimizing the expectation of the mean squared error 
E e (Xf X, t) j with respect to D and S and inter charging 
the order of differentiation and expectation, 

E fxgCz, z, t)] - ])e[x^]-se[xzJ = 0 

and E [Xg (X,i, t) ] - EE [xxj- S E [x^^ 0 

(3»48) 

A stationary random process is orthogonal to its derivative 
14 j that is, 

E [z i 1 = 0 (3.49) 

Hence solving for D and S, 

D = E [Zg (X, Z, t)] / B [F ] 
and S = E [z g (Z, Z, t)]/ E [ Z^ ] (3,50) 

If the nonlinearities in stiffness and damping are linear- 
ly related, g (Z, Z, t) can be expressed as, 

g (z, X, t) = f^ (z) + fg ( 2 ) (3.5l) 

Since, for a Saussian process equation (3.49) ingjlies 
independence of Z and X, from equation (3.50) and (3.5l), 
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and S 


D = E [ X ( Z )]/ E [ F ] 
S = E f X f, ( Z )!/ E r 1 


(3.52) 

(3.53) 


Erom 1iie fl^iores (?) and (s) it is observed that, ia the 
range of relative displaeemeat betweaa axles and mspension, the 
spring chsiacterstics are Iniear, that is, (z) = . Substi~ 

tuting this value in equation ( 3 . 53)5 obtain S = k. To deteimine 
the equivalent linear damping coefficient D equation (3.52) will be 


used. 


given by. 


# 

Eef erring to figure (6), the nonlinear |^8 J f 2 ( X ) is 


fg ( Z ) = a^ Z - bg 

So X 


a . Z + b . 
1 1 


- 00. ^ Z 4 - h2 

- h2 4 i 0 

0 ^Z ^ h^ 

h, < Z < 56 (3.54) 


Equation (3*52) can be written as. 


Z f ( Z ) p ( Z ) d Z 
o 

Z^ p ( X ) d z 


(3.55) 


TOiere p ( Z ) is the probability density hinction. Since 

• f * 

X is a Gaussian xandom variable (zero mean and variance g ; 
p ( i ) is given by 

* o 


(X ) 


xV2 cr^ 


/Vrf (T 


(3.56) 
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Substituting equation (3.54) in (3.55) 

-h 
I 

D 


2 


2 (a^ X - bg) p ( X ) d X + j" Z (g^ X) p ( X ) d Z 

■— b 


h 


1 


^0 


X (g^ z) p ( X ) d X + ^ (a^ X + b^) Z p ( Z ) d I 

(5.57) 




W<r 


Since p ( X ) is an eveia function 

'C;0 


a 


D 


/CO 

l)*2/* * ^ 1 I * * • 

2 y Xp(x)dX+— yXp(x)dZ+ 


r 

h 


(T 


1 


'60 


z p ( Z ) d X + 


^OO 
b 


2 

(T" ^2 


Z^p(x)p(z)dX + 


_2 

2 ^ 

(T ^2 


,*r / ,*, N / * N • ^2 f *2 

Zp(z)p(x)dX +-4 ‘ 


(T ^0 

00 *2 2 

f . - Z V2 r 


.00 

Lot I. = I Xp(x)dX 


Z p ( X ) d X 

( 3 . 58 ) 


h 


/Tf r 


Z e 


d Z 


h 


2 


Putting; Y = 7^/2 (f , 

, 2 /„ 2 
- h /2 (p- 

^ ^ 2 2 

. -XV2£r • 

Z e d Z 


(3.59) 


Let I 


h 


/W r" 


Putting Y = Z / yr (T 


2 / v yr* <r 2 -7“ 


2 (T ( 


Y e dy 
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Let the fb Hewing integral be considered. 


e ^ dy 

9 


If u = e and v = y, then, 


2 / 


e djr = ye'' + I2ye‘' dy 


2 -y 


2 -y" 


or y e " dy 


2 , 2 


Therefore, I, 


2 ^ w/Tfi ' 2 h/ /r<r 2 


jiL ^ 

/tT "" " /w 


Vv^cr 2 


dy) - 


e"^ ay 


hCT' -41^2(1^ 


H erf (^ ) - ^ 


(3.60) 


Similarly, 


* 2 . 

^ ? e-^=<r dX 


r 


2 2 2 
cr ^ f ^ \ ^ b^ -Jti / 2 0 " 

_ — ^ ej;.f0 — ) + : : ® 

2 /Tr m 

Sibstituting, equations b.59), (3.6o) and (3.6 1) m (3.58), 
-h^ / 2 d" ^ 

B = ^ (b^ 4-a^h^ -g^b^) 

-b^/ 2cr^ 

^ ( bg + a^bg - g2 bg) 


vW^T" 


a-, V ^2 ^ / 2 s 

+ — ^ erfc ( — ZTT;) 2 

2 ^ / 2 <r 

g h^ §2 , ^2 s 

+ erf'. ( -HU ) ■*■ 2~ T^T 

2 ynr ^ 


(3.61) 
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Considering the fecial case, when the nonlinear term 
consists ol Coulomb friction and linear viscous dan^jing, that is, 
putting = hg = 0 in equation (3*62), the 83 q)ression for D becomesj 


D 


_J 2 1 °2 


y^TT ^ 


(3.63) 


3.8 DElMffililATIOl'T OP VARIANCE OP VELOCITY 
ACROSS THE mSHPOT 


The variance of the velociiy of sprung mass relative to 
2 

unsprung mass is needed in order to determine the equivalent 

X 

linear damping of the su^ension da^pot. The velocity of sprung mass 
relative to axle is givm by, 


X = Z - W 
Wlieru, Z = Y + L 0 


(3.64) 

(3.65) 


The two terms Y and L 0 on the right hand side of 

equation (3.65) arc random variables with zero mean and variance 

and L® 0"^ respectively. Therefore Z is also a "trend free” 
Y 0 


random variable and variance of Z is given by 27 J 


(TT 

z 


2 2 
= t: 


- (T. 


( 3 . 66 ) 


Y © 

where is correlation coefficient between Y and ©. Similarly 
the mt«n value of X is zero and variance of X is givoi by, 


(f7 


cT. + cn 

z 


-2/ (p. ' (p 
OT Z w 


(3.67) 
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vvhcre ^ ±s the correlation between Z and W » can be calculated 
from ttie following equation, [a] 



R. . ( 0 ) 

w z 


/E/(oy E. (b) 
w z 


( 3 . 68 ) 


vdiere E, . (o) is the cross-correlation Amction between the random 

W Z . . 

variables W and Z. E. (o), E.(o) and E, . (o) can be calculated 

W Z W Z 

from, matrix of cross spectral densities. Applying equation (2,6), 



I 

fOO 



E. (O) = 

1 

i 


( AO) dOO 


Tf 

J 

f w 





-oo 



and E, (o) = 


(03 ) dCO 

(3.69) 

Z 


z 




-00 



Similarly, 


fOO 



R. . (O) 

ss, 


. (03) dOO 

(3.70) 

w z 


/ f 

z 



wiiero S. , (03) is the cross spectral daisity function between 
W Z 

• • 

Y and Z. 



CHAPTER IV 


EESPONSE OP THE VEHICLE TO ROAD ROUGHITESS 
4.1 INIEDIUCTION 

As an automobile runs on a road, it is subject to random 
excitation. The responses, such as vertical displacement and accele- 
ration of different components of the automobile and its passengers, 
are of interest to the automotive engineer. In this chapter, the 
rosponot of the automobile has been investigated. A mathematical 
model of the vehicle is proposed in the last chapter. Por a statio- 
nary road input and a .linear system, the response can be determined 
from the complex frequency response aad the input power spectrum. But 
for a nonlinear system, the response is difficult to determine. How- 
ever, for a weakly nonlinear system, an equi-valent linear system can 
be found out, which in effect should give approximately the same 
reoponse, as would a nonlinear system. Por such a linear system, the 
congjlex frequency response is a function of the input. Response quan- 
tities sought in tills investigation are the displacement of sprung mass 
relative to ttie axles, absolute acceleration of the seats and riders. 

It is also intmded to study the effects of various parameters of the 
system on different responses and to select such a set of parameters 
that would give maximum comfort to riders. The method to determine 
various responses is described In the following sections. 
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4.2 RESPONSE OP A MULTI-DEGREE OP lEEELOM SISTM 


The set ol inputs to 1iie system are represented by a 
vector P (t) and outputs by X (t). The relation between input and 


output is given by, 
t 


X (t) = j ["h (t-zjjF (Z) dT 
h ' 


(4.1) 


where the integration is inteirpreted in the saise of mean square con- 
vorgence and fh (t)] is the matrix of impulse response flinctions, A 
typical elanent h.. (t) in the matrix of impulse response function, 

jK 

describes the motion at node j due to a unit impulse excitation at 
node k applied at time t = 0. The matrix of impulse response func- 
tion and the matrix of complex frequency response form a Pourier 
transform pair ,, 


60 


h (t) = 2^ \ [h ((iJ) J exp (iCOt) doJ 

H (to) * j (-icOt) dt 

^c-0 


(4.2) 

(4.5) 


In tile tiieory of random vibration, the first and second order statis- 
tical properties are of importance. They are determined from the 
relation (4.I). As operators expectation and integration are linear, 


E 


E 


(4.4) 


X (t)] = j [h (t -T)] {Dl dt 

[I ct,)x. p I" [h (t, -zpj.ji erp} 

O 0 

|~h' (tg - 'TTg)"! d d Z" 2 (4-»5) 
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Tho excitation is "trend £ree". Therefore} E j*" X (t) J is 
a null vector, as can be readily seen from equation (4.4), Eor a 
weakly stationary excitation, the coirelation matrix of the excitation 


vector, depends only on the time difference. The matrix is denoted by 


L . P ' ^ 


E 


i (Tg) { 


(4.6) 


If it is assumed that every elemaat in matrix R_ is continuous, ab- 

E 

solutely integrable and of bounded variation in the entire domain of 


the argument, then 


Til cc, -tj)' 

L p 


R (T 1 Iburier transformed; 

up -1 -24 

00 


ft (-)] exp f iCh(t^ -Tg) ^ dCAl (4.7) 

fX) ' ^ 

^ 2^ i Pp ^ ^ J itOu) du (4,8) 

—60 

where (CO ) ^ is tae matrix of cross spectral densities of the 
excitations. Substituting equation (4.7) in equation (4.5) and in- 
terciianging the order of integration, 

X (tp X» (tg)]^ J^H (.-D, t^)] [j2L (o))J [h* (o:>, t2) ^ 

-<Xi " ^ ^ 

exp I i<ft’(t^ - tg)? dO) (4.9) 

-j ,t -j 

H (0>, t) = J jh (u)J exp (~iCOu) du 


E 


whore 




exp (-iOJu) du 


(4.10) 


•CO 

Since the excitation begins at t 


= 0, j^h (u)^ 


is a null matrix for 


t 0. If t and % tend to infinity, keeping the difference t^ - t^ 
constant, equation (4.9) yields correlation matrix of the response. 



tX) 

-JJ- 

exp i CO (t^ - tg) j dCO 
Compfiricic ocjia-bion (4,1l) wi-fc (2,5), 

[t Cw)] = [h (w)][€ (oa)J[H* (u.)] 


(4«1l) 


reduces to, 


(4.12) 


For a single d^ree of freedom system equation (4.12) 


(W ) 

X 


H ((!)) I (CO) 


(4,13) 


■vdiero 


H ( OJ ) is knowi as traasmittancy function. 


4.3 RESPONSE OF TIO-EBSREE OF FREEDOM SYSTEM 

Let X (t) and X* (t) be the displaconent and acceleration in- 
put to the two degree of freedom systoa. The corresponding one sided _ 
power spoctiuffi Sj^ ( DO ) and S^* (CO) are defined ia the range ^ to 0)2- 
For a utationory input the two power spectrum are related: 

S„ ( (J» = t.J' S, ioO) for W, . (2.27) 

X 

Lot Y (t) bo tlic absolute acceleration of the moss II (figure 5). 

Tho frccpcncy response function for X (t) as input and Y (t) as output 


2,2 , 2,2 


are derived earlier. 

i 2 5 2^ ^2^ ^1 ^2 (3,42) 

^ 2I 

where, - i 2 (l +f ) *^2 [^1 


H 


## 

Y 


+ fjL?. 00/ 


2 


( 3 . 41 ) 
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To ij'ot '.iiiput-output power spectrum relationship, 
cqua lion io f-'ppliod : 




H . ( C.t) ) 


S , ( a)) 
X 


(4.14) 


kii tlic input to the linear system is ’’trend free", tiie 


output uli’.o iia!,,; auio moan. The variance can he determined as. 


r t.o 
E ' Y 




/Ocp 


H , (O) ) 
Y 


S,, ( U^) dOJ 
X 


(4.15) 


In order to noud:uuonciDnalize, equation (4.15) is divided by 


E 




H Y 


E I Y J 

E jx^ 1 / ^2 

I, , S., (CO)d(;d 

J(>>1 X 


^ s,, (w) dO) 


(4.16) 


wilCTC. 


oxxd 


i) 


,4.v4 < 2:2 2 .4, .4 

^1 ^2 ^ ^1 ^2 


CO 


(4.17) 


+ (l+J ) U 2 j + LO^ «)/ j ■ 

r -2 to’ (1 ) Sa'^s + 2 1 2 ‘-'^1^ j 


+ 

2 2 


(4,18) 


SirndiciTly, the ratio of variance of relative displacement 
of maso II vdth respect to mass I and variance of displacaient xnput 
can similarly be derived : 





bJ1r_§1-L_ 

B ] 


/ ^2 




Srs I ' («) a« 

A 


S^ ( 6u ) d cO 
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'•i/r' 


(/■ 


H 


RS 


D 


S,, (CV) dO.? 


(4.19) 


iX>2 


a..'. 


S (W) dCO 


#icre H 


^RS 


U' 


( 4 . 20 ) 


4.4 RIDE GOt.ilDRT EVALUAITION 

Tho respcKise of tlie seats inside ttie vdaicle acts as the 
tojwt to the passenger. The input vibration is described in teims of 
pow(.r opectam. Eospmse of human heii« to random rtbration Input has 
not boun ahidied mudi. Howsvar in the low freguenoy range response of 
persons to sinusoidal inputs has been extensively studied. Eesilts 
obtained by various investigators differ from each other. Ihis nay be 
due Id variations in phyelologieal reactions In Indiviauals as well as 
fatigue (phyolool os well as mental) level of test eubjeote at the time 
Of oonduoting, the experiment. In this investigation, the recommendation 
offered by the International Standard Organisation (I.S.O., 
dered. (figure 12) 

Eocently it has be«. investigated Aether die evaluation 

1 . - Tar Q-iniisoidal vibration, could also be 

system applicable only to sntpl 

• -j. P T- ml-xtures of sinusoid-al vibrations, 

used with well fbunded validity for 


T2^mln| 



"'”1 

1 

i 

. 4 J 


! ! 

i 

i 

, 


limit multiply acC^j 
leration values by 2 

.1 i « * 


For reduced comfort 
divide acceleration 
values jby 3.15 


4 6 1 10 20 

Frequency. - Hz 


40 60 80 


too 


0 PROTOSEO : RECOMMENDATIONS FOR THE 
“ECTS OF VERTICAL VIBRATION ON MAN 
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and for rmdoiii vibrations imposed on seated persons [28 | . With 
rofcroncc to tolerance 1 imit it has been observed that formation of 
a general effective ■'/alue from the square root of the sum of squares 
of tile weighted accolcjration amplitudes of two vibration components 
does not exactly correspond to the subjective physiological effect. 

In this investigation the man is assumed to be a linear 
syf?teu, heviiy; r, transfer function defined as the ratio of ride sensa- 
tion to input spectrum. The transfer fiinction T (f) is a function of 
frequency f, find is given by 

T (f)-lMl} 

^oro, 2Z (f) is ride sensation 

Z (f) io input to the rider. 

Lot C (f) be the numerical value of the r.m.s. acceleration, 
a huni£si'i beint-; caii tolerate (figure 12, 8 hour limit curve). The 
tKUi;:*i-er lUnction T (f) can be determined from the following relation, 

0 (4' 


T (f) 


C (f 


(4.22) 


where, G (4) is tlie niuiuffic^l value of r.m.s, acceleration at frequency 
of 4 Hz. Erom equations (4.21 ) and (4.22) equivalent nde sensation 

at of 4 IIz cm 0. deterntoed. If 8 (f) - 

tlUJi to ttio rider, tte output usrlanoe can be determmed from the 

follov/Jng integral, 


'-2 

'1 


Var * \ & (f) 


(4.23) 
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4.5 OP EESPONSES BY mraiGA-I 

INTBORATION OP SYSTEM BIPPEElEIifTIAl EQJATIOHS 

In Section 4»2, a method is described in whicii the responses 
c..re obviined uoing frequency response function. Instead of non- 
linear pf-xnaetors, equivalent linear system parameters are used in 
tile calculations. The validity of calculation in frequency domain 
analysis cf.n be checked by an alternate approach of determining 
the response of tiie nonlinear system by direct numerical int^aration 
of tiic system differential equations. 

The system is described by 6 second order differential equa- 
tions. Prom, these second order equations, 12 first order differen- 
tial equations are obtained using the following notations : 

Y^ (t) = (t) 

Yg (t) = Eg 

(t) - y (t) 

y^ (t) = 0 (t) 

y^ (t) = (t) 

\ (^) » Zg (t) 

(t) = (t) 

\ = K 

Yg (t) = Y (t) 

Y^Q('t) = ® 

y^^(t) = z^(t) 

and Y^gCt) « Zg (t) 
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liquation (4.24) in (3.14), (3.15) and (3.18) to ( 3 . 21 )^ 

T, (t) 


^2 

(t) 
Y. (t) 
Y 5 (t) 

i (t) 


(t) 


(t) 
Yg (t) 

^ 1 ,( 1 ) 


M 


1 

3i 

M. 


Y^ It) 


^9 


■1-^1 -""l 

Yg (t) + (t) - Y^ (t) ] 

\j2-\ ^2 


CD] 


X 

*2 


M. 


[Yj (t) + \ (t) - Yg (t)] 


i [yj + Ii I4 - q ] -r [h* S - ^ih 

n 

-[ 


X 

M 


[ 


q + q q 


q * q qo - q 


K, 

T 

K, 

“I 


-q] -1 

[q^qq-q] 

a[Y,.i,q-q] tq^qqo-'X 


Y - Y,^ 
10 11 


] 

1 


qo‘^’ 


®. I*. 
1 1 


®2 ^2 


K. 1 . 

1 1 


"^2 ^2 


q ■" ^1 q ■ j 

hb 

~ j 

_q ■" ^1 ^10 ■ q J 

q + q q - q_ 

\ ' ^ ■ 
d j 1* «i r 

q + q qo - q ] 

^3 "*’ ^1 ^4- ” 

wo» 

1 1 

”■ J 

fq + q qo - qq 

— * 

T ^2 ^2 ] 

j “ J ; 

jjg + Ig - Y^g] 
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(t) 


Js - ^^3 

” 1. 

Y, ' 

S : 

Y 

~ Y - 1 

Y 


1 

1 

4 

“^1 i 

_ 11 

9 1 

10 ^ 

Ct) 

= 

- \ 

- lo 

Y. ' 


Y „ 

~ Y„ - 1 

Y 

liip 

_o 3 

2 

4 

lUg 

12 

9 2 

10 


(4.25) 


Tlioyc; cquationc can be numerically integraied, 
uocd ifourtl: or>ler laetliod is 1iie me credited to Gill 


The most widely 



'3 



["3.1 " = - 


J 

yr 



vihere 


^ 3,3 




fj (ic, ) 


(4.26) 


f. 


} m • I 


K 


3,4 




• « », 


yp” 

l/ C, 


h K + (l + ^<o -z I 

^ ^S2,2 ^ ^ 12,3 J 


Using (4.26), equation (4.25) can be dntegrated. The values 

of Y^, and Y. are stored ia each iterations. Subsequently the 
3 3 3 

variance of Y,, 1, Y. are detemiaed and compared wi1h those ob- 
tatned by onploying equivalent linearization technique. 
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4.6 mi’A OP A Ji'ASSaiCrER CAE 


Tht: xrilbTrifVbion regarding a passenger car, supplied by 
jriiuiuutiui r.o Icrr;; Ltd. |8 ] ere given below : 


Total 


Pront 


Hear 


2563 lbs. 1401 lbs. 1162 lbs, 
3043 lbs. 1538 lbs, 1505 lbs. 


3363 lbs, 1619 lbs. 1742 lbs, 

148 lbs. 276 lbs, 

1390 lbs, 1229 lbs. 


1, a) Curb weigl.t or unladen 

b) Lorui'i la dun 

(l Puitnungor in front, 

2 Per.amguro in roar, 

+ 3 y. 20 Ibo. luggage) 

c) Rilly laden 
(bCfO lbs, mxiiiiuiQ) 

2, Unupiun^^ niass 

3, Spnuig ma.uiJ 

4, Monuait of Inertia of vebicle 2 

pitching about its C.G, is 1900 slug ft 

5, Wlitiul bast; 8‘ 1" 

6, Gtjat location in vehicle : 

Rouraoat 16 5/8" awsiy from Eear Axle centre. 

Pront tjont 48 13/16" away from Eear Axle centre. 

Ihc ouapension spring and damper characteristics and static 
lc«d dcilection curve fbr tires are shown in figures (7) to (ll), 

4.7 EEaJLTS 

The response of two degrees of freedom system is investi- 
gated using B^pations (4.16) and (4.19). lixe tire stiffness can not 
be varied arbitrarily* weight of unsprung mass at 
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dou!! not oiTi.T' ituch choice for variation of them. Uaially the rear 
fixlu io henvioi' tiitu the front, because it houses the differential 
etc. Tiio p*'r! meters that can be varied are suspension spring stiff- 
nec;;, d'-aiphjg value and prtion of spiung laass tiiat each subsystem 
hue; 'to carry. The variable parameters are^V^ (natural frequaicy of 
os.-concl mss), 5*2 coefficient of second mass) and f (ratio 

ol' i.iunril and maos l). The ratio of variance of acceleration of mass 
IJ to thi'l of acceleration input for varying values ofCt^, f ^ 

dfctcrminod using equ? tion (4.16), These values are plotted in figu- 
res 13, 14 and 15. Similarly, the latio of variance of displacement 
of bwohU to displacement input variance are found out by means of 
equation (4.19). These values are plotted in figures 16, 17 and 18, 


The equivalent linear damping of the dashpot with nonlinear 
cliartictcristics are found out using equation (3.62). in iteration 
pTOQCwa ia used to achieve convergence. The input to the vehicle xs 
khoviai in too fora of power spectrum. The road rou^ness can be varie 
by ciumgtog toe expression for power spectrum. This can be achieved by 
varying too rougtoess coefficient (see ^2.1.5.). The variation of 
equivalmt linear damping with respect to roughness coefficient is 
determined and plotted in figure (l9). further, the responses are 
dependent on toe equivalent linear damping. The responses corres 
ponding to different road conditions are plotted m figure(2o). 


, data toe roughness coefficient can be 

Sbr a hnovaa road uara, me 

-pnT the svstem is determined 

detorminod. Ooa^lex frequmey respo 
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unin/; c;(|intion (3.44). The response of the vehicle for a given 
road condition can be found out with the help of equation (4.12), 

Ihoi't; arc several paiometer 'vdiich can he varied to change the res- 
ponse quantities. The effect of variation of suspension spring 
stiffness, suspension damping, seat spring stiffness, seat damping 
and location of centre of gravity of sprung mass, on response of 
vdiicle ere given in Tables 1, 2, 3, 4 and 5 respectively. 

Transfer function of human being as a function of frequency 
is calculated from the I.S.O, proposed recommoadation for the effect 
of vertical vibration on Man (figure 12, 8 hrs curve), using equation 
(4.22), The power spectrum of the rider response is the product of 
transfer function of the rider and power spectrum of the seat response* 
Variance of rider acceleration (corresponding to 4 Hz frequoicy) can 
be obtained from Uie power spectrum of rider response. By varying 
■tile suspension parameters, ■ftie variance of seat and rider acceleration 
•is "varied. The variance of seats and rider acceleration are shovan. in 
Table 6 for various seat and damping -values. 

-6 

Ibr a road vd'bh roughness coefficient 32 x 10 , random 

data is generated. The method ■tn generate data from power spectrum 
is discussed in article 2,4. These data are used as input to tiie 
system. The systan differential equation (4.25) is numerically 
tegrated by Gill's modification of Runge Kutta method. The algori-thm 
for -the integration technique is given by equation (4.26). The 
■variance of acceleration ofbpuncc motion is computed and is found to 
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TABLE 1 

EEIECT OE VilRIATION OE SJSPEUSlOEf SPRIITa STIEEEJESS 01 

TEHIOLE EESPOISES 

= 1200 0^ = 1400 =60 dg = 60 = 100 Eg " 

Front sus- jRear suspsn-jYariance of iVariance of jYardance of ivariance of 
pension spr-|sion spring Ibounce acclnfpitcii accln^jifront seat xrear seat 
ing stiffne^stiffness } p 4 I -4 o 2 , 4 

Its, /ft, I lbs, /ft, 1 ft, /oec. I Sec. 1 ft. /sec, a ft /sec»_ 


1100 

1100 

10.52 

.087 

3.91 

3.98 

1100 

1300 

11.36 

.082 

4.60 

4.76 

1100 

1500 

11.60 

.084 

4.91 

5.29 

1300 

1100 

11.24 

.087 

4.39 

4.25 

1300 

1300 

11.66 

.083 

5.04 

5.11 

1500 

1100 

11.31 

.094 

4.49 

4.38 

1500 

1300 

11.54 

.091 

4.86 

4.86 

1500 

1500 

12.23 

.087 

5.35 

5.76 
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mSLE 2 

EFFECT OF VARIATION OF SUSPENSION Ili'ffl>ING ON VEHICLE RESPONSES 


= 1300 Sg = 1500 0^ = 1200 = 1400 = 60 dg = 60 


Front sus— J[Eear sus— p/arian oe of Xvarisnce of 5 Variance of 
pension dam-Jpension dam-|bounce acclnfpitch accln.ffront seat 
ping -value |pjng value I o a I _a faccln. . 

Ib./ft./sec.H lb,/ft«/sec][ ft. /sec, I Sec, f ft«^/sec. 


j Variance of 
Xrear seat 

laccln.^ 

( ft « /sec, 


4 


125 

100 

11.87 

.089 

5.45 

6.01 

125 

125 

12.84 

.075 

5.23 

5.44 

125 

150 

13.61 

.074 

5.10 

5.23 

150 

125 

13.63 

.073 

5.21 

5.35 

150 

150 

14.46 

.067 

5.13 

5.16 

150 

175 

15.26 

.068 

5.14 

5.20 

175 

150 

15-31 

.065 

5.25 

5.13 

175 

175 

16.16 

.061 

5.30 

5.18 

175 

200 

16,97 

.063 

5.41 

5.37 

200 

175 

17.03 

.059 

5-51 

5.18 

200 

200 

17.88 

.057 

5.64 

5.38 

200 

225 

18.68 

.058 

5.81 

5.67 
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nsm 3 

EI'i'ECT OP YmU.TLO'E OP SERIEG STIPPtTESS 01 THE 7PHIGIE 

EESPOISE 

= 1300 Sg = 1500 = 100 Dg = 120 =60 dg = 60 

Praat seat |Rear seat | Variance of ^Variance of j[ Variance of Xvariance of 

spring stiff Ispring fbounce acdafpitch accln.ffront seat |rear seat 

-ness fstiffness I o 4 ^ 4 laccln.g ^|acGln. ^ 4 

lbs. /ft, \ lbs. ft, \ ftf/sec, \ Sec. \ ft. /sec.L_. 


1100 

1100 

1 1 -167 

.086 

5.42 

5.53 

1100 

1200 

11.92 

.082 

5.47 

5.62 

1100 

1300 

11.95 

.082 

5.46 

5.67 

1200 

1200 

11.94 

.082 

5.54 

5.65 

1200 

1300 

11.96 

.082 

5.53 

5.67 

1200 

1400 

11.99 

.083 

5.53 

5.71 

1300 

1300 

11.98 

.082 

5.60 

5.67 

1300 

1400 

12.00 

.082 

5.60 

5.71 

1300 

1500 

12.02 

.083 

5.59 

5.75 

1400 

1400 

12.01 

,082 

5.67 

5.71 

1400 

1500 

12.03 

.082 

5.66 

5.75 

1400 

1600 

12.05 

.083 

5.65 

5.80 
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table 4 

EFFECT OF VARI/ 1 .TIOF OF D/.MPEIG IH SEfi.T OR VESIOBE HBSPORSE 
= 1300 Sg = 1500 C^ = 1200 = 1400 = 100 = 120 


Front seat ] 
damping | 

Ib/ft./sec. ^ 

[Rear seat fVariance of ][ Variance of 

[dairping {bounce accln][pitch accln, 

[ 1 J 

i lb,/ft,/sec{ ft^/sec.'^ { Sec.~^ 

{variance of {i^riance of 
, {front seat {rear seat 

4 ““^ 2/ + 

1 ft. /sec.i ft, /sec. 

40 

40 


12.13 

.092 

5.86 

6.21 

60 

60 


11.99 

.082 

5.55 

5.53 

80 

80 


11.68 

.078 

5.32 

5,42 

100 

100 


11.41 

.074 

5.24 

5.26 

120 

120 


11.16 

.071 

5.25 

5.20 

140 

140 


10.94 

.069 

5.30 

5.20 

160 

160 


10.73 

.067 

5.39 

5.23 

180 

180 


10.54 

.064 

5.48 

5.28 

200 

200 


10.37 

.062 

5.59 

5.36 

220 

220 


10.21 

.061 

5.71 

5.44 

240 

240 


10.06 

.060 

5.82 

5.53 
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TABLE 5 

EEFE3T OP LOO ATI Oil OP OIHTRB OP GRAVITY OP SPEUKG IltASS 
OK THE VEHICLE RESPOKSB 



s^ = 1300 

s^ = 1500 

= 1200 

1 

= 60 

' °2 = 
"2 = 

1400 

60 

= 100 

Dg = 120 

Listance iDistnace jListance o: 
of front ][of rear ][ front seat 
axle from|axle frani|froni C.G. 
O.G, of xC.G, of xof sprung 
sprung 5 sprung Imass 

mass ft. ]imass ft. ][ ft. 

fiDistance jYariance i^YariancejiVarianceiVaraance 
|of rear fof touncelof pitch][of front? of rear 
fseat fromfaccln, ](accln. Jseat gseat 

|c.G. of I 1 liaccin. |accln. 

fsprung Ip 45 - 4 I ,45 „^2 , 4 

tmass ftlft /sec. I Sec. ? ft./seci ft^s^ 

3 

-5 

-1 

-3.67 

11.77 

.126 

5.27 

5.53 

3.25 

-4.75 

-.75 

-3.42 

12.02 

.102 

5.49 

5.64 

3.5 

-4.5 

-.5 

-3.17 

12,00 

.087 

5.53 

5.69 

3.75 

-4.25 

-.25 

-2.92 

11.97 

.080 

5.52 

5.73 

4.0 

-4.0 

.0 

-2.67 

11.91 

.078 

5.48 

5.73 

4.25 

-3.75 

.25 

-2,42 

11.84 

CO 

0 

♦ 

5.41 

5.70 

4.5 

-3.5 

.5 

-2.17 

11.76 

.095 

5.32 

5.62 

4.75 

-3.25 

.75 

-1.92 

11.66 

.112 

5.21 

5.50 

5.0 

-3 

1.0 

-1.67 

11.54 

.136 

5.10 

5.30 
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TiffilE 6 

STUDY OP Y/iRIATIOI OP VARIANCE OP ACCEIERATION OP RIDER WITH 
RESPECT TO VARI/iNCE OP ACCEIiERi TION OP SE/-T 

= 1300 Sg = 1500 C^ = 1200 Cg = 1400 D^ = lOO D 2 = 120 


Pront and iVariance of 

rear suspen- ffront seat 
Sion damping iacoln* . 

16, /ft, /sec. I ft /sec.^ 

1 Variance of 
irear seat 
](accln.„ . 

1 ft /sec 

1 Variance of 
Jaccelera toon 
lof rider on 
Ifront seat- 
I , ftVsec.'^ 

^Variance of 
Jacc deration 
|of rider on 
|r ear seat i, 
t ftvsec,^ 

40 

5.87 

6.21 

3.65 

3.86 

60 

5.53 

5.71 

3.43 

3.53 

80 

5.32 

5.42 

3.33 

3.38 

100 

5.24 

5.23 

3.35 

3.33 

120 

5.25 

5.20 

3.38 

3.35 

140 

5.30 

5.20 

3.46 

3.39 

160 

5.39 

5.23 

3.55 

3.46 

180 

5.49 

5.28 

3.64 

3.54 


200 


5.59 


5.36 


3.73 


3.63 
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2/4 

be 25 *89 ft /sec. Using equiT/alent linearization technique j ihe 
•variance of acceleration of bounce is found to be 24.96 ft^/sect 
which compares well with that obtained by using numerical integra'td-on 
technique, 

4*8 SELECTIOl OF OPTIl-IOM PiBAHSEMS 

In the data given on article 4.6, -the moment of inertia of 

sprang mass about its centre of gravi'ty is not Icnown, It is assumed 
2 

to be 763 slug ft . This satisfies -fche dynamic decoupling condition 
(refer equation 3»25). After decoupling into two subsystems the 
following quantities for a passenger car (§ 4,6) are calculated* 


OABIE 7 

COldPUTED CHARACTERISTICS OE A PASSENGER OiE 

FRONT REA-E 


Mass I (m^) 

2.3 slugs 

4.28 

slugs 

Mass II (m^) (noimal laden) 

23.9 slugs 

23.8 

slugs 

Natural frequency of Mass I (60^) 

70.5 rad/sec 

51.7 

rad/ se ( 

Natural frequency of Mass II (^ 2 ^ 

7.55 rad/sec 

8.4 

rad/se< 

Ratio of Mass II to Mass 1 {'f ) 

10.4 

5.55 


Ratio of natural frequencies of 
Mass II to that of Mass I 
( CO 2 /CO^ 

.107 

.162 


If 'the vehicle moves on a road at 60 m.p.h. wi'th rough- 
ness coefficient 32 x 10"^ (equation 2.17), the damping values for 
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f!ront and rear wheels are found to be 250 and 201 Ib/ft/sec respec- 
tively (using figure 19)* The corresponding damping coefficients 
can be calculated j using equation 3*339 and ere found to be •? and 
• 51 respectively* 

Prom the figure 14 ^ the ratio of variance of acceleration 

of mass II to variance of acceleration of disturbance input ( (T^^/ (p^) 

y s; 

—4 —4 

are found to be #56 z 10 and #52 x 10 « The frequency range of 

interest for human vibration is 1 to 80 cps* Integrating the power 

spectrum of displacement road input (using equation 2.2?) in this fre- 

queacy range the variance of input acceleration ( — •^*) is found to be 

5 2 , 4 

2.885 X 10 ft /sec. . Therefore, the standard deviation of accelera- 
tion of mass II (<r^*) of front and rear portions are 4 and 3.86 

y 

respectively. 

Similarly the standard deviation of relative displacement 
of mass II with respect to Mass I ( SZ) a3?e found (refer figure 17) 

Hu 

to be ,208" and .224". If the clearance be taken as 5 SVo #iidi 

turns out to be 1,04" and 1.12" respectively, Ithe probability that the 

displacement: of mass II relative to mass I will be within these -values 

is 0.9999997. But in actual vehicle the clearance pro-vided for -the 

front and rear suspensions are 4" and 3.75" in normal laden condition 

and 3.5" and 3" in fully laden condition (figure 7 and 8) vdiich are 

far above even -the safe 5 (fT„ -vaiues. Thus the clearance seems -fco 

ES 


be overdesigned 
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It can be seen that for the actual bump stop distance provi- 
ded ^ the vehicle can ihoto on a more rough road at the speed of 60 m#p*h* 

-6 — ' 

e#g* an average principal road with roughness coefficieat of 128 x 10 

It is rx)t desirable that the spring should be compressed to 
its solid length. The 5 value is taken to be 2.5’^* The variance 

of displacement of mass II relative to mass I is given by^ 

^TeS "" ^ '5 x^'l2 ^ .00174 ft . 

The variance of displacement for the road roughness can be determined 
by integrating power spectrum for displacement road input (using equa- 
tions 2.6 and 2,17) in ihe frequency range of 1 to 80 Hz. This is 

found to be ,00187 ft^. Therefore, 

_ 2 

<ris ^ ,00174 g. 

^2 .00187 " ' 

^ X 

Ibr the values, ,6, ~^ = 10,4 and 5.5 for the front and rear 

portions, the corresponding iA }2 / (jJ^ found from figure 17. These 
■values are .0875 and .15, correspondiag to C*)^ ■values of 70.5 and 51.7 
(table 7), the -values of are 6.17 and 6,75 rad/sec respectively. 

Using equation (3,53) the front aid rear suspension stiffness are 
calculated to be 920 and 1090 Ibs/ft. 

Many of "the data required such as seat spring and dashpot 
characteristics are not a-vailable. If "these da^fca are hnown, the 
optimum parameters for maximum comfort of rider can be determined. 

The effect of "varia-tion of different parameters on the system per- 
formance will be discussed in the next chapter. 



CHAPTER V 


GONCnj SIGHS 

5.1 CEHEEAL OBSERYi'TIOlT AFD nn mTmT 'TS 

J^oia results obtained in the previous chapters j following 

observations and coniments can be ciade ; 

(i) The acceleration of Mass II (pig. 5) increases as Odg / 
increases (Pig, 13, 14). Keeping (jJij constanty acceleration 
of mass II increases as OO 2 increases. Prom equation (3.33)? 
it can be observed that acceleration increases as stLffiiess 
increases. Since mass II does not offer mch flexibility for 
a change? the decrease in acceleration can be achieved by 
decreasing the stiffliess. But the decrease in stiffness 
causes increase in the displacement of mass II relative to 
mass I (Pigi. 16? 17 and 18). Since the bump stop distance 

is fixed? the stiffness of suspension can not be decreased 
arbitrarily. If the spring is coirpressed to its solid length 
the spring acts as a rigid link between axles and chassis. 
Thus the rider is subject to a very hi^ value of accelera- 
tion causing discomfort. 

(ii) Acceleration of mass II increases as danping coefficient ^ 
increases (Pig, 15). Thus for better comfort it is re- 
commended to have as low damping as peimitted. But damping 
can not be lowered below the optimum road holding values fl6j' 
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(iii) Figure 15 shows iiiat as ^ (ratio of mass II to mass l) 
decreases, the acceleration of mass II increases. There is 
little scope fcr varying the unsprung mass I but mass II can 
be varied to a certain exteat bjr placing the centre of gra- 
vity of sprung mass so that equal acceleration is obtained 
at both front and rear portions of the vehicle. However, 
the position of O.&. between the two axles can not be varied 
to a great extent and thus it appears that the parameter ^ 
does not offer much scope for variation. 

(iv) The suspension danping characterf tics are nonlinear in nature. 
The equivalent linear damping value is determined applying 
equivalent linearization technique (equation 3*62). From 
figure (19) it can be observed that the damping value de- 
creases as the road roughness coefficient is increased. 
Equation 3,63 confirms this result. Physically it can be 
explained as follows. The relative velocity between Sprung 
and unsprung mass increases as the road roughness Increases, 
The damping in shock absorber consists of linear viscous 
damping and Coulomb friction. Although the linear viscous 
damping fcrce increases proportionately with relative velo- 
city, Coulomb friction does not increase correspondingly. 

Thus the overall damping force does not increase propor- 
tionately to road roughness. At high value of road rou^- 
ness, the friction force in dashpot is very small compared 

to the linear viscous damping force. Therefore, the damping 
value remains almost constant. 
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As the road roughness is doubled, the acceleration 
response of sprung mass and seats does not increase two fold 
(Pig. 20). This is due to decrease in damping force. How- 
ever at higher -values of road roughness it increases propor- 
tionately since the damping value remain almost constant. 

(v) As expected, if the spring stiffness of front and rear 
suspension is increased, the bounce mode increases (Table l), 
but the pitch motion is not much affected. 

(vi) As damping in suspension increases, the bounce loode increases 
(Table 2) because as damping increases more of disturbance is 
transmitted, 

(vii) AS seat spring stiffness increases, bounce acceleration and 
acceleration of seats increases (Table 3)» It is desirable 
to employ a spring with a low stiffness. The criterion coula 
be fixed by a given s-tatic deflection viien the rider sits 

on it. 

(viii) As seat damping is increased from a -very lov/ value, accelera- 
tion of sprunjj mass toad seat decreases (Olable 4). If ttie 
seat dampjxig is fartb.er increased beyond a certain value 
the seat acceleratioii in creases • This is because for low 
values of damping, an increase in dampxog value reduces the 
oscillations of seats where as for high values of damping, 
an increase in damping helps in transmitting more of /ibra 


tion from spruaog mass to seat* 
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(ix) If the centre of gravity of sprung mass is shifted ffom 
centre of axles, the bounce and seat acceleration may 
decrease slightly as can be seen, from table 5, but pitch 
acceleration increases. Because pitch motion causes more 
discomfort than bounce, the centro of gra.vity should be as 
near as possible to the centre of two axles. An increase 
in vheel base will cause decrease in pitch oscillation. 

(x) Ihe equivalent acceleration of rider at 4 Hz is calculated 
according to equetion (4.2l). Prom table 6, it can be seai 
that the variance of acceleration of rider varies as variance 
of acceleration of seat. Therefore the former need not be 
calculated except vhere it is desired to design the seat 
spring and to know the comfort conditions. 

5.2 AWUES POH POETHPE RESEARCH 

Ptor off road conditions, the assumption regarding stationariiy 
of ground profile cay not be justified. In such cases, -ttie concept of 
continuously changing power spectrum or "evolutionary power spectrum" 
oan be used. Expressions for tiie response quantities for the vehicle, 
can be developed in tarns of the evolat ionary ^ectrum of the terrain. 

There are various sources of inputs to the passaigei and 
driver. Por instance vibrations are transmitted to the driver through 
steering besides seat aid foot rest. If ihe vibration transmitted 
through the steering oolumn can be minimized, the handling character- 
sties of the car can be considerably ijaproved even for rough roads. 
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To improve the ride comfort, study can be cade of the vdiicle with 
seat and foot rest separated from the chassis by spring and dashpots. 

Planes that land on aircraft carriers get a short runv®y 
disiance. They are subject to ver3’' high decelerations. Therefore, 
the responses momentarily spurts up to "very high "v^alue due to decele- 
ration of the vehicle. This case needs further study. 

In this investigation, the effect of engine vibration has 
been neglected. Engine transmits periodic disturbing forces to the 
chasis. At certain speeds, aada as idling, the frequency of disturb- 
ing forces due to engine may coincide with resonant fre^ency of 
vehicle. In this case the disturbing forces due to engine, however 
small it may be, affect the comfort of passengers. The disturbing 
forces due to wind gust and other aerodynamic forces may be investigated, . 

Effect of vibration on man has been studied for determinastic 
inputs, but random vibration teclmiques has not been properly atteap- 
ted. An extensive study of human comfort ^ould be done and mental 
and physical fatigue criteria need be established more realistically. 

Response of large vehicles carrying vibration sensitive caxgo 
or human beings and having core than 4 wheels should be studied. This 
extends to buses, large trucks and ambulances etc. 

Off road vehicles, such as tanks, military vehicles etc., 
need be properly analysed for rmgh terrain profiles. The criteria 
in these cases might not be comfort, hit of proper tolerances, ability 
to align guns to the target or for radar detection etc. 
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